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We consider the consequences for circulating-current states of a cyclic, four-spin, "ring-exchange" 
interaction of the type shown recently to be significant in cuprate systems. The real-space Hartree- 
Fock approach is used to establish the existence of charge-current and spin-current phases in a 
generalized Hubbard model for the Cu02 planes in cuprates. We augment the Hartree-Fock ap- 
proximation to consider correlated states renormalized by Gutzwiller projection factors, which allow 
us to gauge the qualitative effects of projection to no double occupancy in the limit of strong on-site 
Coulomb repulsion. We find that charge flux states may be competitive in the underdoped and 
optimally doped regime in cuprate systems, whereas spin flux states are suppressed in the strongly 
correlated regime. We then include the ring-exchange interaction and demonstrate its effect on 
current-carrying states both at and away from half-filling. 
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I. INTRODUCTION 

The possibility of phases supporting circulating cur- 
rents in two-dimensional (2D) systems has been debated 
for some time. The orbital antiferromagnet, or charge 
flux (CF) phase, was first identified in the context of exci- 
tonic semiconductors"'^ and reintroduced very early in the 
analysis of high-temperature superconductors,^''^ while 
the spin flux (SF) phase was introduced within a general 
planar model for the latter.'* On a square lattice, the CF 
phase [Fig. 1(a)], has counter-rotating currents of car- 
riers on neighboring plaquettes, whence the alternative 
terminology "staggered flux" phase, and preserves SU(2) 
spin symmetry but breaks time-reversal invariance. The 
SF phase [Fig. 1(b)], which can be said to be a fermionic 
realization of a spin nematic, may be considered as super- 
imposed counter-rotating current patterns in which up- 
and down-spin currents are equal and opposite (no net 
charge transport), and preserves time-reversal invariance 
but violates SU(2) spin symmetry. 

A considerable body of theoretical work now exists 
concerning these states, particularly the CF phase. Early 
analyses employed the Heisenberg-Hubbard model with 
an SU(A^) generalization of the spin in order to obtain a 
controled mean-field decoupling, ^■'^ while further studies 
have since been performed of the different saddle points 
obtainable in the t-J model. ^ The relationship of the CF 
phase to the superconducting one has been explored in 
detail within an SU(2)-symmetric slave-boson gauge the- 
ory of the the t-J model by Lee and coworkers,^ and 
by using a Hubbard AT-operator formulation.'' More gen- 
eral models for the cuprate planes, including not only 
a finite on-site Hubbard interaction (J7), but also ex- 
plicit nearest- neighbor superexchange (J) and Coulomb 
interaction [V) terms, have been considered using weak- 
coupling techniques. At half-filling and for [/ > 0, 
the Hartree-Fock approximation* (HFA) and more so- 
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FIG. 1: Schematic representation of (a) charge-flux (CF), and 
(b) spin-flux (SF) phases. Solid arrows represent currents of 
spin-up electrons and open arrows spin-down electrons. 



phisticated renormalization-group analyses of divergent 
susceptibilities®'^ agree on a phase diagram with mag- 
netic order for weak J and V , a charge-density-wave 
(CDW) phase for strong F > 0, CF for F > and suffi- 
ciently large J > 0, and SF for F > and J < 0. Nega- 
tive V (nearest-neighbor attraction) yields superconduct- 
ing phases of singlet symmetry if J > and triplet sym- 
metry if J < 0.*" 

Interest in the CF phase has been revived recently by 
an accumulation of experimental evidence suggesting the 
presence of a quantum critical point in the superconduct- 
ing region of the phase diagram.** One of the propos- 
als for the "hidden" order parameter associated with the 
transition is the conventional CF phase, also known as 
the d-density-wave (DDW) state, *^ although more com- 
plex current patterns based on a three-band description 
have also been proposed.*^ To date the only experimental 
means of identifying a CF state is the breaking of time- 
reversal symmetry, and the interpretation of one such 
measurement using polarized angle-resolved photocmis- 
sion spectroscopy*'* remains open. Circulating-current 
phases have been considered numerically in the context 
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of exact diagonalization (ED) studies of chiral fluctua- 
tions in the 2D t-J model/^ and fluctuating staggered 
currents on a doped ladder/^ while a more general model 
displaying steady currents was considered by density- 
matrix renormalization-group (DMRG) studies in a lad- 
der geometry. 

The possible significance of cyclic four-spin, or "ring- 
exchange", interactions in cuprates was also recognized 
at an early stage. The important term arises at fourth 
order in a strong-coupling (small t/U) perturbative ex- 
pansion of the single-band Hubbard model, and has 
been shown in this limit to give the leading correc- 
tion to the nearest-neighbor Heisenberg model. Exper- 
imentally, inelastic neutron scattering measurements of 
the spin-wave spectrum in the 2D system La2Cu04,^'' 
and in the quasi-one-dimcnsional spin-ladder compound 
La6Ca8Cu2404i,^^ as well as two-magnon Raman scat- 
tering measurements on doped systems related to the 
latter, indicate the presence of ring-exchange contribu- 
tions. The magnitude of this interaction has been quanti- 
fied as being 20% of the nearest-neighbor superexchange 
interaction, i.e. K ~ J/5, both by systematic expan- 
sion of the three-band Hubbard model, ^'^ and by a self- 
consistent spin-wave analysis yielding good quantitative 
agreement with the experimental data.^^ 

At the qualitative level, the presence of a ring-exchange 
term has been found by consideration of exactly soluble 
extended Heisenberg models to generate a rich variety of 
novel phases on a ID spin ladder. A number of these 
phases have been found and characterized by ED and 
DMRG studies of the isotropic ladder.^® However, for 2D 
systems accurate results beyond the Hartree-Fock anal- 
ysis of static magnetic phases^® remain scarce, and the 
only numerical study to date was performed for a system 
with XY spin interactions.^^ 

The ring-exchange interaction on a square lattice may 
be considered as the product of a virtual hopping process 
around an elementary plaquette, while the CF and SF 
phases are those in which real currents circulate on the 
plaqucttcs. However, to our knowledge there have been 
rather few studies of the mutual effects of fiux states and 
ring exchange. The issue has been addressed for the CF 
phase of a t-J model with unprojected nearest- and next- 
neighbor hopping by applying a large- approach. This 
study suggested that a ring-exchange interaction acts to 
suppress the homogeneous CF phase, but may favor in- 
homogeneous current patterns and other states involving 
a further breaking of translational symmetry.^^ 

Here we adopt a qualitative approach to gain fur- 
ther physical insight into both the real-space nature of 
circulating-current phases on the 2D square lattice and 
the infiuence upon them of ring-exchange interactions. In 
Sec. H we present the generalized Hubbard Hamiltonian 
supporting flux phases in the weak-coupling limit and 
discuss its treatment in the HFA. Sec. HI contains a sum- 
mary of the properties of the CF and SF phases obtained 
in HFA with unprojected hopping. We turn in Sec. IV 
to a study of electronic correlation effects, which are in- 



cluded within the Gutzwiller approximation through pro- 
jected hopping terms, and this allows us to make contact 
with the physical parameter regime for cuprate systems. 
In Sec. V we include the ring-exchange interaction and 
demonstrate its effects on both CF and SF phases at dif- 
ferent fillings by analyzing the plaquette current and the 
ring-exchange energy contribution. Sec. VI contains a 
summary and conclusions. 



II. MODEL AND FORMALISM 

A. Extended Hubbard model 

To establish the existence of CF and SF phases (Fig. 1) 
in a 2D lattice we consider the extended Hubbard model 
in the form 



H. = Ht + Hu + Hj 
where the first two terms 



Hv + H, 



K, 



(1) 



Ht = -t 5^(4c,>+H.c.)-<' ^ (4c,-,+H.c.), 

(ijV ((«.i>>o- 
Hu = U'^niiUii, (2) 

i 

constitute the conventional Hubbard Hamiltonian for 

a single, nondegcncratc band with nearest- and next- 
neighbor hopping. Here 4 denotes the electron creation 
operator, riia — cl^Cia, ni = ^„nia, and the notatation 
(((«, j))) indicates a pair of (next-)nearest-neighbor 
sites. We adopt t = 1 as the unit of energy. The terms 

Hj + Hv = JY,^i-^j + V niarija', (3) 

correspond, respectively, to superexchange and Coulomb 
interaction between nearest neighbors. The t-U-J-V 
model (1) has been considered by a number of authors, 
including a Hartree-Fock treatment of the infinite system 
with only nearest-neighbor hopping^ and RG studies of 
both the half-filled case^ with t' = {] and the model with 
small t' near the van Hove fillings.^ 

Finally, the ring-exchange term in Eq. (1) takes the 
form 

Hk = kJ2 [(S^-S,)(Sfe-SO + (SrS,)(S,-Sfc) 



-(Si-Sfe)(Sj-S;) 



(4) 



where □ denotes the sum over the four spins (ijkl) of 
each plaquette labeled clockwise. The phase diagram 
of a 2D system in the presence of this term has been 
considered only for the Heisenberg antiferromagnet,^^ al- 
though further progress has been possible for ID mag- 
netic systems. ^^'^^ 
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B. Hartree-Fock Approximation 

The HFA is the simplest and the most efficient ap- 
proach for the study of interacting electron systems. Al- 
though its application may seem inherently limited to 
the regime of weak interactions, it provides satisfactory 
results at all interaction strengths when the correlation 
energy is small (for example in magnetically polarized 
phases) , and thus can be used to obtain meaningful qual- 
itative insight into the nature of possible ordered phases 
also for large U. In the HFA, multiparticle interac- 
tion terms are decoupled into all possible operator pairs, 
all but one of which are replaced by c-number expec- 
tation values to leave a single-particle Hamiltonian. In 
the present analysis the choice of operator pairs with fi- 
nite expectation values is as follows. Because our pri- 
mary interest is the study of flux phases, we will work 
with a repulsive intersite Coulomb interaction > 0, 
in which case no pairing decouplings need be considered, 
(cjj^cj^) = = (c|^cj^). In the absence of single-site spin- 
flip terms in the Hamiltonian [Eq. (1)], there is no rea- 
son to expect finite averages of the form {cl^Cj^}, where 
a denotes the spin projection opposite to a. Thus flux 
phases will compete or coexist with magnetic and pos- 
sibly charge-ordered phases, and the finite expectation 
values to be considered are the site charge 



the site magnetic moment 

mi = (Si) = 1 i4a^iocl3Cif3), 
a,/3=T,i 

and the bond order parameters 

Sijcr = (cJo-Cjo-)- 



(5) 



(6) 



(7) 



In a circulating-current state, the latter parameters are 
complex, s 



I] a 



with the real part corre- 
sponding to the bond kinetic energy and the imaginary 
part to the bond current. 

For later reference we present the Hartree-Fock decou- 
pling expressions for all of the interaction terms in Eq. (1) 
in the form in which they will be used here. The Hubbard 
term is approximated by 

Hu-UY^ [{ni^;)nii +ni^{nii) - {nii){nii) 

i 

- i^i CiT ) 4t - (4t ) 4i CiT + (c|^ Ci4 ) (c|^ CiT ) J ,(8) 

where the signs follow from fermionic statistics. We de- 
fine the effective bond order parameter 



9ijcr — ^ijcr ~t~ 2^ij^^ 



(9) 



in terms of which a decomposition of the J, V, and K 
interactions gives the following one-body terms in the 



Hamiltonian matrix. For the superexchange interaction 

+\J Y ct^c,-, + XaS^cla,) (10) 

+ lJ E [{St-iXaS^)clcj^ + {S^-iXaS])clci^] , 

where A| = 1, = —1, and {gija, djia-} are complex 
quantities. For the nearest-neighbor Coulomb interaction 

Hv - VY + rijcl^Cia^ 

— V y ] (^^jio-cl^Cja- + Sija-C^j^Cia^ . (11) 

Finally, the ring-exchange interaction is decoupled into 
single-particle terms multiplied by all possible combina- 
tions of the averages obtained by contracting the remain- 
ing six operators, 



{v).o- {kl){{ij)) 

+ hK J2 [{S^-iXaS])clc,^ + X.S]clc„ 



X J2 (Sfe-S;), 

{kl){{ij)) 



(12) 



with the intersite spin correlation function, determined 
within HFA, given by 



(Sfe-Si) = {Sk){Si) — ^{skiisiki + Skiisiki) 

-idswrl^ + lsfeul^)- 



(13) 



The structure of the decomposition is similar to that of 

the superexchange term, although in this case the net 
contributions are quartic in (Si) and (c|^Cjct'), so may 
be expected to be small at the Hartree-Fock level. Note 
that in Eqs. (10-12) we have omitted the constant terms 
bilinear in the expectation values [analogous to the final 
term in each line of Eq. (8)] . 



C. Self-consistent calculations in real space 

We have performed calculations on a finite cluster in 
real space to establish the microscopic conditions under 
which the CF and SF phases are manifest as circulat- 
ing currents of real electrons (or holes), which arise from 
a mismatch in hopping expectation values between the 
two directions along each bond in a coherent pattern. 
These bond expectation values, the site charges, and the 
site magnetic moments arc computed for any set of inter- 
action parameters and iterated to self-consistency. The 
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HFA in real space has been used extensively for qualita- 
tive modeling of cuprate systems, and in particular for 
the analysis of charge-inhomogeneous stripe and other 
density-wave phases, which it tends to favor. A full re- 
view of the technique is presented in Ref. 29, and we 
summarize here those aspects of its application which 
are important for the current analysis. 

Most importantly, the HFA suppresses fluctuation 
terms and thus favors (homogeneous or inhomogeneous) 
ordered states which minimize static interactions. Mag- 
netic and charge-density interactions, which are zeroth- 
order in particle-hopping fluctuations, are enhanced rel- 
ative to kinetic terms, and it is for this reason that 
the method tends to find ground states with inhomo- 
geneous charge distributions and/or magnetic order. By 
contrast, methods such as the Hubbard X-operator and 
other large-A'' approaches have a tendency to favor ki- 
netic terms, as a result of which they do not show stripe 
phases and often produce flux-phase ground states. Thus 
a real-space Hartrec-Fock analysis is not an appropriate 
framework for an unbiased discussion of energetic con- 
tributions and true ground states; in specific terms, the 
flux phases we investigate are frequently local rather than 
global energy minima. Mindful of this fact, we will re- 
strict our considerations to qualitative statements con- 
cerning the nature of circulating-current phases and the 
relative effect upon them both of doping and of the dif- 
ferent interactions, particularly the ring-exchange term. 

We note also that the method is especially sensitive to 
the choice of starting state, as well as to boundary condi- 
tions and in some cases to conimc^nsurability effects due 
to the electron filling. A careful specification of these ex- 
trinsic variables is required in order to make meaningful 
comparisons yielding intrinsic electronic properties.^'' We 
will consider a system of size N = 12x12 sites with pe- 
riodic boundary conditions; qualitatively similar results 
were obtained with systems of size 8x8 and 16x16, and 
we have established that flux states may be found for 
any rectangular system with even-length sides. Open 
boundary conditions lead to a suppression of circulat- 
ing currents and to inhomogeneous edge values, and are 
not considered further. The establishment of global, as 
opposed to local, energy minima may be pursued as in 
Ref. 29, but for comparative studies of flux states is not 
of prime importance (preceding paragraph). Finally, by 
considering the evolution of the results with electron fill- 
ing we will demonstrate that commensuration effects are 
absent and thus that the system size is sufficient for ro- 
bust qualitative conclusions. 



III. CIRCULATING-CURRENT STATES 
A. Dependence on doping 

We begin by establishing homogeneous flux phases in 

the generalized Hubbard model on a cluster. As expected 
from the weak-coupling results mentioned in Sec. I,^'*'^ 
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FIG. 2: Doping-dependcnco of kinetic energy T and current 
/ per bond for each spin orientation as obtained in HFA for 
CF (J > 0) and SF (J < 0) phases with U=\J\=AV = 4t. 



we find that the nearest-neighbor interactions J and V 
are of crucial importance: circulating-current states ap- 
pear close to half-filling only for simultaneously large val- 
ues (> 1) of the interaction ratios \ J\/U and AV/U . The 
CF phase, with finite currents of up- and down-spin elec- 
trons moving in the same direction [Fig. 1(a)], occurs 
when J > and SF states, with counterflowing up- and 
down-spin currents, appear for J < 0. The parameter 
sets for flux-phase formation are dictated by the strong 
propensity within the real-space Hartree-Fock approach 
towards charge localization, driven by large values of U , 
charge inhomogeneity driven by large V ^ and phase sep- 
aration driven by large J. In this section we characterize 
the dependence of the homogeneous CF and SF phases 
on the doping and on the magnitude of the interactions. 

We define the filling as the normalized number of elec- 
trons in the system, n = Ne/N, and the hole doping as 
X = \ ~n. We choose the parameters U =\J\ = AV = At, 
and show in Fig. 2 the doping dependence of the bond 
order parameter [Eq. (7)] for a single spin species and a 
single direction along an x- or a y-axis bond. The real 
part T = Re Sijc is a component of the bond kinetic en- 
ergy, and constitutes one eighth of the hopping kinetic 
energy per site of the system, while the imaginary part 
I = Im Siju corresponds to one half of the bond current 
for a single spin species (and thus to one quarter of the 
charge current in a CF phase). It is clear that, for mod- 
erate interaction strengths, circulating-current states are 
robust close to half-filling. However, the current shows 
a monotonic decrease as a function of doping until the 
CF and SF phases become unstable at doping values of 
0.16-0.18. 

The magnitudes of the (negative) real parts Tqf and 
TsF are scarcely affected by the presence of the currents. 
In the absence of flux phases or other types of order, T is 
expected to be constant at low doping for the cluster un- 
der consideration, beginning to drop beyond x = 0.153; 
this number does not affect the behavior of established 
ordered phases. For the intermediate coupling strengths 
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FIG. 3: Properties of CF and SF phases with increasing in- 
teractions U =\J\= 4:V for doping x = 0.028 as obtained in 
HFA: (a) kinetic energy T and current / per bond for each 
spin orientation; (b) total energy E and kinetic energy ST. 



and the second line of Eq. (11). Figure 3(b) compares the 
kinetic energy contribntion to the total energy, illustrat- 
ing the dominance of the static, zeroth-order interaction 
terms contained in Eqs. (10) and (11). While the kinetic 
contribution remains nearly constant, the total energy 
increases monotonically, reflecting the penalty for fail- 
ure to adopt charge- or magnetically ordered states. It 
is this failure which ultimately destabilizes the CF and 
SF phases, rendering them uncompetitive in comparative 
HFA studies in real space. 

One observes again that growth of the current has lit- 
tle effect on TcF in Fig. 3(a), and that Icf saturates with 
increasing interaction strength, approaching the limit 
I^cfI = |-^cf| where the phase per bond is Xija = 7i'/4. 
Thus the phase (j) for a plaquette approaches (f) = tt, and 
the CF state is the 7r-flux phase of Ref. 2. By contrast, 
in the SF state Isf continues to rise with increasing in- 
teraction strength, the kinetic contribution Tgp falls, and 
this phase becomes unstable to a symmetry-broken mag- 
netic state with finite net moment. We note also that the 
flux states remain stable in the limit of vanishing interac- 
tions, although in this case they are not the global energy 
minima of the system. A similar situation is found by 
analysis of different saddle points in a large- approach: 
for the S\]{N) saddle point, which yields the CF state, 
the circulating-current phase also survives at zero inter- 
action strength. While Fig. 3 illustrates the situation for 
X = 0.028 (4 holes in a 12x 12 system), the same behavior 
is found at higher doping: for x ~ 0.1 the evolution of 
the bond phase in the SF state is almost identical, while 
for the CF state it converges more slowly to 7r/4, and the 
current in the zero-interaction limit is somewhat smaller. 



used in Fig. 2, the kinetic contributions in the CF and 
SF phases are suppressed when the currents are maxi- 
mal. This reduction is by approximately 5% below the 
value found in the order-free "metallic" phase realized 
beyond the critical doping x ^ 0.17. We comment on 
the fact that the current "order parameter" , Iqf , drops 
continuously to zero with increasing doping, while IgF 
changes discontinuously. This contrasting behavior may 
be taken as a reflection of the different symmetries bro- 
ken by the two phases (Sec. I), through the coexistence 
with competing phases allowed by these symmetries. 



B. Dependence on interaction strength 

The CF and SF states for small doping are character- 
ized in Fig. 3(a), which shows the quantities T and / 
as functions of interaction strength for a fixed choice of 
the interaction ratios \J\/U = AV/U = 1. Qualitatively, 
the ciirrent does not contribute to the overall kinetic en- 
ergy, and the finite imaginary part of the bond hopping 
parameter is stabilized by Fock terms originating from 
decouplings of the J and V interactions, as seen in the 
presence of the negative signs in the first line of Eq. (10) 



C. Dependence on next-neighbor hopping 

The effect of a next-neighbor, cross-plaquette hopping 
t' in Eq. (2) is shown in Fig. 4. We have computed 
the bond-order parameters for both nearest- and next- 
nearest-neighbor bonds for charge- and spin-flux phases 
over the range of values of t' for which these phases are 
maintained. The qualitative behavior is the same for all 
dopings in Fig. 2. and is illustrated for an intermediate 
choice of x. We have deflned t and t' with inclusion of a 
negative sign in Eq. (2), so t = 1 and t' <Q corresponds 
to the physical sign for cuprates. We note also that for 
clarity of presentation Fig. 4 shows |r| — ~T but the true 
sign of T' . The most interesting feature is that small val- 
ues of t' have no effect at all on CF and SF states, in that 
the expectation values of the bond energy and current, 
which are determined largely by Hj and Hy in Eq. (1), 
are not changed by the additional kinetic term. In this 
regime, which persists more strongly for negative t' , di- 
agonal hopping gives a very small contribution to the 
energy which is positive for t' < Q. Negative values of 
t' show a small tendency to reinforce both types of flux 
phase, by which is meant to increase the bond current 
with a concomitant small decrease in bond energy, which 
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\t'\ < 0.3, whereas the SF phase is entirely suppressed. 
Further studies of the role of diagonal hopping, including 
consideration of circulating-current phases established on 
the diagonal bonds at very large values oi\t'/t\, have been 
conducted in Ref. 28. 

We conclude this survey of the qualitative properties 
of flux phases in the weak-coupling regime by comment- 
ing that in the unrestricted Hartree-Fock approach es- 
sentially no inhomogeneous current- carrying states were 
found. While the method presents no barriers to phase 
coexistence, over most of the parameter range the co- 
existing flux-phase components arc small compared to 
other, competing phases, and only a very weak, uni- 
form current pattern is maintained. This situation may 
be altered in the presence of a relatively large next- 
neighbor hopping parameter t', which also raises the 
possibility of obtaining specific inhomogeneous current 
configurations,^^ but this parameter regime is of little 
physical interest in the context of cuprates, and is not 
considered here. 



IV. FLUX PHASES WITH STRONG 
CORRELATIONS 



FIG. 4: Influence of next-neighbor hopping t' on flux phases 
with U = \J\ = AV = 4t for doping x = 0.083 as obtained in 
HFA: kinetic energy T and current I for each spin orientation 
per X- and j/-bond, and kinetic energy T' per (x+y)- and {x- 
j/)-bond, are shown for (a) CF and (b) SF states. 



in the context of cuprate systems is a noteworthy result. 
Positive values of t' show a more pronounced tendency to 
suppress both types of flux phase; for the CF state there 
is a stable small-current regime over a range of values of 
t', in which the next-neighbor hopping contribution to 
the kinetic energy is quite significant. 

The regions of stability are cut off for both signs of t' 
as the need to minimize the total kinetic energy drives 
transitions to phases with similar energetic contributions 
from both nearest-neighbor and diagonal bonds. These 
states are marked by a sharp rise in \T'\. For the CF 
state this process occurs for \t'\ ~ t, and involves homo- 
geneous zero- or small-current phases. For the SF state 
the transitions involve inhomogeneous phases of small to- 
tal kinetic energy and non-zero magnetic moments, and 
occur for significantly smaller \t'\. This latter result may 
be taken as evidence that the SF phase is generally less 
stable to competing phases [Fig. 3(b)]. We stress that the 
values of t' required to destabilize the CF and SF states at 
this filling should not be considered as quantitatively rep- 
resentative, and emphasize only the qualitative behavior 
in the presence of t' . The results are consistent with the 
renormalization-group analysis of Ref. 9 for the physical 
sign of t' [t' < in our notation convention) , in which the 
location of the CF phase boundary is scarcely affected for 



A. Gutzwiller Approximation 

The analysis of the preceding section is limited to 
the regime of weak coupling. However, interactions in 
cuprates are strong, and in particular the large on-site 
repulsion U causes a suppression of double site occupa- 
tion due to electronic correlations, which is not taken 
into account in the imrestricted HFA. As a consequence, 
the parameter regime in which fiux phases were found 
in Sec. HI has little relation to the model parameters 
appropriate for superconducting cuprates. In this sec- 
tion we use the Gutzwiller Approximation (GA), which 
includes strong correlations through their eflPect in sup- 
pressing states including doubly occupied sites from the 
Hilbert space, and thus allows us to approach the cuprate 
parameter regime. Because the GA reproduces qualita- 
tively the effects of (partially or fully) projected hopping 
on the single-particle level, an adapted HFA of the type 
employed in Sec. Ill remains valid. Recent considera- 
tions of Gutzwiller wave functions for multiband models 
have also led to an improved understanding of the Mott 
transition in these systems.^" 

While many of the early analyses of flux phases in 
cuprate-related models employed an eff'ective large-f de- 
scription such as the t-J model, relatively few of these 
employed the projection of the hopping term associated 
with the constraint of suppressed double site occupancy. 
This constraint, which with full projection yields a hop- 
ping term of the form 
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t' V [{l-ni^)clcj,{l-nj^)+R.c.], (14) 



E 



is technically very difficult to implement systematically. 
Because of its singular nature the connection between 
unprojected (independent electron) and projected models 
remains quite unclear, and the possibility arises that flux 
phases may indeed be far more relevant for the latter. 
One exact means of implementing projected hopping is 
in numerical studies, such as the exact diagonalization 
analysis of uniform and staggered flux phases in the t- 
J modeP^ and of the CF phase in the presence of hole 
pairing.^^'^^ 

Here we approximate the projected model by the 
Gutzwiller approach, ''^ in which electron correlation ef- 
fects are contained within multiplicative prefactors mod- 
ifying the terms of the Hamiltonian. These prefac- 
tors are obtained from statistical considerations based 
on the numbers of states available and excluded by the 
constraint. This approximation to the effects of a large 
U term was presented for the 2D Hubbard model in 
Rcf. 36, where the technical details are also summarized. 
For the purposes of controling the extent of projection, 
here we apply these factors in their general form at finite 
U (see also Ref. 37) 
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2d{n - 2d) 



^22 



md^ 
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where gt, gj, and g™"* multiply the corresponding 
terms in Eq. (1). The variable n is the electron filling 
(Sec. HI. A), and the parameter d represents the density 
of doubly occupied sites in the system, d = {ni^riii). Be- 
cause we consider only those uniform phases found in 
the HFA analysis of Sec. HI, the density d is indepen- 
dent of i. For finite values of U, three terms are required 
to represent the nearest-neighbor Coulomb interaction, 
depending on the occupation state of the sites involved. 
Because the number of bonds with at least one doubly 
and one singly occupied site is low for all parameters and 
dopings under consideration, it is a very small and thus 
acceptable approximation to proceed with a single factor 
gv = gv- 

The effective model with projected hopping is consid- 
ered to be derived from a Hubbard model with large 
but finite U.^^ Thus the parameter d is nonzero, in ac- 
cord with the finite values of the J and V terms, whose 
physical origin lies in virtual hopping processes involv- 
ing an effective double site occupancy. We define the 
quantity do (n) as the average value of the double site oc- 
cupancy for a fixed, finite value of U at each electron 
filling n, as determined in the HFA to the correlated 
state {do{n) < (n/2)^ because this reference Hartrec- 
Fock state is magnetic). An effective t-J-V Hamiltonian 



with projection factors gt, gj, and gv determined using 
do(n) in Eq. (15) should then represent the starting t- 
U Hamiltonian. Within this framework, projection to a 
lower net double occupancy d < do{n) is a variationally 
motivated procedure which seeks to minimize the total 
energy by further elimination of the energetically costly 
doubly occupied sites, the limit of this process being a 
full projection with d = 0. In Sec. IV.C we will use d as 
a variable parameter to illustrate the interpolation be- 
tween the unrestricted HFA and the strongly correlated 
state for the same electron density n. We note, however, 
that only small corrections to the energy and magnetiza- 
tion are expected due to this procedure in the strongly 
correlated regime, because the unrestricted HFA already 
captures most of the relevant effects, and for large U ap- 
proaches the same limit as the projected model. While 
this partial projection procedure is of limited meaning in 
the context of a model with arbitrarily chosen values of 
J and V as in Sec. Ill, the effective model used in this 
section is that obtained by large- [/ expansion of the Hub- 
bard model (2). The parameters J and V in Eq. (1) thus 
have the specific values J = At^/U and V = — J/4. 

Further insight into the physical meaning of the GA 
and its relation to the Gutzwiller projection for cuprate 
systems, may be obtained from the consideration of vari- 
ational wave functions, a study which has been under- 
taken both for superconducting states^^ and for coex- 
isting superconducting and CF states. "^^ We note here 
that the partially projected states used in this framework 
are not appropriate for addressing the Mott transition at 
half-filling, as has been shown systematically from their 
violation of optical sum rules. At finite doping these 
considerations bear a certain similarity to the recent dis- 
cussion of "gossamer" superconducting states. ^^■''^ 

For the purposes of analyzing flux phases in a cuprate- 
relevant model with projected hopping, we begin by con- 
sidering a Hubbard model with U = 12. The Hartree- 
Fock density do (n) of doubly occupied sites is then com- 
puted for all values of the filling 0.8 <n<l(0<a;< 
0.2). An effective projected model is obtained by com- 
puting the Gutzwiller factors [Eq. (15)] at each filling by 
varying d between the values d = 0, which describes full 
projection and is implemented to extract the properties 
of correlated electrons as described by a t-J model, and 
d = rfo(n), which should reproduce qualitatively the un- 
projected (Hartree-Fock) states. The final version of the 
effective model is then as in Eq. (1), but with t replaced 
by gtt, J replaced by gjt/3, V by — gyi/12, and no on- 
site repulsive term. We emphasize that the derivation of 
this effective model involves a singular truncation rather 
than a systematic transformation, with the consequence 
that no meaningful energetic comparisons are possible 
between the results of Sec. HI and the results obtained 
after projection. However, it is important to understand 
the general properties of the projected system, and com- 
parisons performed within it remain well defined. 
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B. Qualitative consequences of Gutzwiller 
Approximation 

Qualitatively, one of the principal effects of projec- 
tion is to sharpen the competition between the different 
candidate ground states. The exclusion of doubly occu- 
pied states eliminates most of the possibilities for a co- 
existence of localized and delocalized phases, and phase 
boundaries would thus be expected to become better de- 
fined. For CF phases these boundaries may also move 
closer to a parameter regime meaningful for cuprates (see 
below): when viewed as a means for dynamical holes to 
avoid each other (by circulating around each other on 
plaquettes), the CF state indeed allows a greater spatial 
separation than the superconducting state while retain- 
ing a favorable kinetic energy. However, application of 
this schematic picture for the SF phase, which requires 
in principle that up- and down-spin particles can move 
past each other on each individual plaquette, suggests 
that this type of state would be completely suppressed 
by projection. 

By inspection of the Gutzwiller prefactors (15) it is 
clear that close to half-filling the properties of the pro- 
jected model are significantly different from those of the 
unprojected model. With full projection the factor gt 
vanishes as x ^ 0, giving a complete suppression of 
the kinetic contributions to the total energy, and also 
of any bond currents. This suppression is large for small 
doping, where the kinetic and current terms are signifi- 
cantly smaller than in the unprojected system. However, 
because all competing states arc; similarly affected by 
projection, this does not imply that circulating-current 
states must become less competitive. By contrast, for 
small d the prefactor gj 4 as x 0, suggesting a 
relatively strong enhancement of magnetic interactions. 
Once again, because the superexchange interaction has 
both spin and bond-order decouplings [see Eq. (10)], this 
trend does not necessarily imply a suppression of flux 
states. The prefactor gv is smaller than gj, and is not 
expected to play a significant role. For larger values of x 
these statements remain generally true while U is large, 
and d correspondingly small, but the relative strengths 
of the kinetic and magnetic terms will change almost lin- 
early with X. 

C. Partial projection 

We begin the more quantitative part of the analysis 
by considering the evolution of the flux phases with the 

extent of projection of the hopping term. Following the 
above discussion we define a partial projection variable 

p = rf/do(n), (16) 

such that p = corresponds to fully projected hopping 
(14), and p = 1 to a projection at finite U which rep- 
resents (by reference to the average number of doubly 
occupied sites) the HFA of Sec. HI. Fig. 5 shows the 
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FIG. 5: Kinetic energy T and current I per bond for each spin 
orientation in a CF state for with partially projected hopping 
for X = 0.028 and J = -W = t/3. 



bond energy and current of a CF state function 
of p for a system with 140 electrons (x = 0.028); while 
this doping illustrates most clearly the effects of projec- 
tion, qualitatively identical results are found for higher 
values. At larger values of p, which correspond to re- 
laxation of the constraint and thus to a larger relative 
kinetic contribution, ITcf] and Jcf change linearly with 
p, and have a fixed ratio (corresponding to a plaquette 
phase (/) ^ 2tt/3) which decreases with increasing doping. 
As full projection {p = 0) is approached the suppression 
becomes stronger, and the real and imaginary parts tend 
to fixed values determined by the filling; their ratio in- 
creases with projection up to an equivalent flux </> = O.8871" 
for this doping. 

The SF phase exists only for negative values of the su- 
perexchange J, which is not a valid parameter regime in 
the framework of a small (t/U) expansion for the Hub- 
bard model. If the sign of J is simply inverted, the cal- 
culation does not converge to a stable local minimum 
corresponding to a SF state for any partial projection. 
The projection factors contribute in two ways to this 
result: the fact that 7sf approaches ITsf] {(j) t^) re- 
duces a difference term [which for the CF state is a sum, 
causing the difference visible in Fig. 3(b)] contributing a 
negative energy, and the suppression of the total kinetic 
energy removes what is then an essential negative contri- 
bution. This failure to form an SF state in the projected 
real-space HFA, which extends to all values of the doping 
X, may be taken as evidence that the schematic picture 
of counter-flowing up- and down-spin currents, which im- 
plies that electrons must hop past each other on each site, 
is in fact a valid description of this phase. 

With regard to the total energy of partially projected 
CF states, the increasing kinetic contribution as projec- 
tion is relaxed {p ^ 1) means that the fully projected 
state has the highest energy. We have commented in 
Sees. II and IV. A on the restrictions concerning meaning- 
ful energetic comparisons because of both the nature of 
the HFA and the singular projection scheme. Within the 
manifold of projected CF states, changes in total energy 
as a function of p are dominated by the kinetic energy 
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FIG. 6: (a) Doping-dependencc of kinetic energy T and cur- 
rent / per bond for each spin orientation in the fully projected 
(p = 0) CF phase with J = -4V = t/3. (b) Corresponding 
bond-order parameters Re Sij and Im Sij . 



term, and, in the framework of a consistent derivation 
from the Htibbard model, no systematic study of ener- 
gies, for example as a function of interaction strength 
[Fig. 3(b)], is then possible. However, it is clear from 
our calculations in the projected regime that the energy 
difference between the CF state and other competing 
phases, such as antiferromagnetic order, is small. 



While we caution as before that quantitative con- 
clusions cannot be drawn, and note also that energy 
scales are altered by projection, it does appear mean- 
ingful to state that, unlike in the unprojected frame- 
work, CF phases are strongly competitive in a projected 
model. This statement applies not only to the most fa- 
vorable choice of interactions (Sec. Ill) but specifically 
for cuprate parameters. The enhanced competitiveness 
of the CF state indicates that projection provides a better 
reflection of the energy gained by dynamically stabilized 
phases, the origin of which lies primarily in bond-order 
decoiiplings of the J and V terms. Wc may conclude 
from this consideration of the GA that CF phases are 
not in fact unrealistic for cuprate interactions when ap- 
propriate projection of the hopping terms is performed. 
More quantitative analyses of this issue are required. 



D. Dependence on doping 

We consider next the dependence of the fully projected 
{p = 0) CF state on doping. The evolution of the bond 
parameters |Tcf| and Iqf with x, shown in Fig. 6(a), 
exhibits three important qualitative features. First, the 
effective projected model approaches the tt-Aux phase, 
|Tcf| = IcF, when the doping decreases towards half- 
filling. However, in this limit the kinetic energy and cur- 
rent also vanish in the projected model, in contrast to 
the unprojected one. Secondly, the kinetic component 
increases linearly in magnitude (becomes more negative) 
with hole doping across the full range of x, tending to- 
wards the value found for a projected metallic state (not 
shown). Finally, the bond current grows at first, but be- 
yond X ~ 0.08 falls continuously to zero (at Xc — 0.16), 
indicating the instability of the CF state to a metallic 
(Fermi-liquid) phase at higher doping. 

We stress that in the projected model the kinetic terms 
are no longer identical to the real and imaginary parts of 
the bond-order parameter Sij , but are related by a ratio 
of Gutzwiller factors (15). The components of s^ , which 
are determined by the bond decoupling of the terms Hj 
and Hy in Eq. (1), are shown as functions of doping 
in Fig. 6(b). The absolute value of the bond amplitude 
\sij\ decreases only slightly with doping, while the bond 
phase falls continuously from 7r/4 at half- filling to zero at 
Xc- At a qualitative level, the monotonically decreasing 
Im Sij and dome-shaped form of Iqp are rather simi- 
lar to the behavior of the pseudogap amplitude and or- 
der parameter for rf-wave superconductivity computed in 
the projected BCS wavefunction framework in Ref. 39. 
Given that we have chosen to allow only bond order, 
but no pairing decouplings, such a similarity may be ex- 
pected from the SU(2) symmetry relationship of CF and 
superconducting phases in this class of models.^ The mi- 
croscopic origin of a coexistence of paired and CF states 
has been investigated by a number of authors. At 
the quantitative level, the net hopping kinetic energy per 
site \Et\ — SITcfI in Fig. 6(a) is for all dopings within 
15% of the value calculated in Fig. 4(a) of Ref. 40 us- 
ing a projected wave function approach to the SC phase. 
This result underlines the energetic similarity of the SC 
and CF phases, and the need for a detailed and unbiased 
analysis to distinguish between them. 

We discuss only briefly the dependence of the CF phase 
on the strength of electronic interactions in the extended 
Hubbard model. Within the framework we have estab- 
lished for the Gutzwiller analysis this means an alteration 
in the value of U, and for increasing U one approaches 
the asymptotic regime where do{n) oc (t/U)^. However, 
as do{n) becomes smaller than its value for ~ 8, the 
Gutzwiller projection factors (15) are governed primarily 
by the doping x. This reflects the fact that increasing 
doping causes a gradual release of the constraints aris- 
ing from strong correlations. In the large-f/ regime of 
most interest it is thus the filling rather than the dou- 
ble occupancy which is the most important parameter 
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FIG. 7: Influence of next-neighbor hopping t' on the projected 
CF phase with doping x = 0.083, showing kinetic energy Tgf 
and current 7cf per spin for each x- and y-bond, and kinetic 
energy T^p per [x+y)- and {x-y)-hond. 



determining the current and the bond phase of the CF 
state. 



E. Dependence on next-neighbor hopping 

Following the analysis of Sec. Ill, we consider also the 
effect of a next-neighbor hopping t' on the fully projected 
CF phase, shown in Fig. 7. As observed also for the 
unprojected model (Fig. 4), small values of t' have no 
effect on the bond-order parameters T and /, and the 
T' contribution, which is positive at t' < 0, is essentially 
a residual consequence of the prevailing CF state. The 
regime of stability for the projected CF phase extends 
considerably further for negative than for positive t' , with 
its midpoint at t' ^ —0.3: it appears reasonable to state 
that for the physical magnitude and sign of t' in cuprates, 
next-neighbor hopping does not disrupt and may in fact 
reinforce the CF state. 

The CF phase is again terminated by abrupt transi- 
tions to zero-current states for sufficiently large values 
of t'. For t' > a transition occurs at rather small t' 
to a phase with strongly reduced bond kinetic energy. 
However, beyond this we find a further CF phase with 
a small but finite current, which permits a significant T' 
contribution. This new phase persists to t' /t > 1 before 
being replaced by another state, of zero current, which 
again has strongly reduced \T\ due to manifestly frus- 
trated hopping. For t' < the conventional CF phase 
persists down to t'/t ^ — 1 before being replaced by a 
zero-current state which achieves negative kinetic energy 
on both nearest-neighbor and diagonal bonds,^^ with \T\ 
only slightly suppressed from its value in the CF state. 



V. RING-EXCHANGE INTERACTION 

We turn next to the effects of the four-spin ring- 
exchange interaction on the flux states of Sees. Ill and IV. 
Because the structure of this term is quadratic in bilinear 
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FIG. 8: Effects of a ring-exchange interaction K on unpro- 
jected CF and SF phases for U = \J\ = AV = At at doping 
x = 0.083: (a) current I per bond per spin orientation; (b) 
energy contribution Ek- 



spin terms (4), the structure of its Hartree-Fock decom- 
position (12) is similar to that of the superexchangc term 
(10). However, in addition to nearest-neighbor spin and 
bond expectation values it contains also next-neighbor 
bond-order terms corresponding to the diagonals of each 
plaquette. These appear in the Hamiltonian matrix in 
the same position as next-neighbor (t') hopping terms. 
The quartic nature of the ring-exchange interaction K, 
and the relatively small expectation values of the order 
parameters ((Sj-S^) < 0.25, |sjj| < 0.25), suggest that 
the effect of this term in an unprojected model will be 
small for all but the largest values of K. 



A. Unprojected model 

The qualitative influence of the ring-exchange interac- 
tion on the improjected CF and SF phases of Sec. Ill may 
be assessed by considering the changes in bond currents, 
and also the total energy contribution, due to the pres- 
ence of the four-spin coupling term. Figure 8(a) shows 
the bond current in a system with a representative in- 
termediate value of doping for a canonical choice of flux- 
phase interaction ratios U =\J\= AV = At over a range 
of values — | J| < K < \J\. It is clear that a positive ring- 
exchange interaction suppresses the CF state, whereas it 
acts to enhance the current in the SF state. 

These qualitative results may be explained by consid- 
ering the "interference" of currents on opposite bonds of 
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FIG. 9: Doping-dependence of Ek in the presence of a ring- 
exchange interaction K = \J\/5: (a) for unprojected CF and 
SF phases with U = \ J\ = W = U; (b) for the projected CF 
phase. 




FIG. 10: Effects of a ring-exchange interaction K on the pro- 
jected CF phase with J — — 4V = t/3 and doping x — 0.083: 
(a) current / per bond per spin orientation; (b) energy con- 
tribution Ek,cf- 



each plaqucttc, which arc coupled by the K term in the 
manner shown in the first hne of Eq. (12), and specifi- 
cally from the role in this interference of currents with 
opposing spins [second term of Eq. (13)]. The suppres- 
sion of the CF state by a ring-exchange interaction is in 
agreement with the findings of Ref. 28; enhancement of 
the SF state is also consistent with a recent exact diag- 
onalization study of a small cluster with K as the dom- 
inant interaction, which revealed a type of SF phase. 
We comment further on this result below. 

The ring-exchange contribution to the total energy, 
shown in Fig. 8(b), is by its quartic nature always pos- 
itive {oTC K > 0. This energy, Ek, constitutes a signif- 
icantly smaller penalty in the SF state than in the CF 
state. Examination of the total energy (not shown) re- 
veals a similar increase for the CF state, reinforcing the 
result of the preceding paragraph. For the SF state the 
total energy also shows a small increase of the same or- 
der as Ek, suggesting that changes in the kinetic and 
magnetic energies induced by the enhanced current do 
not compensate the positive ring-exchange contribution. 
As a function of filling, the influence of the K term on 
the quantities characterizing the flux phases remains as 
shown in Fig. 8 and displays no strong variation over the 
range < x < 0.2. 

We have also computed the bond currents and ring- 
exchange energies of flux phases with K = \J\/5 for the 
full range of doping values over which these exist. The 
current differs only very slightly from its values at if = 



(Fig. 2), although there is a slight suppression which 
brings the CF transition to a lower doping value. The 
energy Ek, which is illustrated in Fig. 9(a), shows only 
minor changes with increasing x until the transition to 
a different (non-flux) state brings abrupt alterations at 
the critical dopings. The relative changes in current and 
energy arc both understood readily from the evolution of 
Sij in the flux states at if = (Sec. III). However, we 
draw attention of the magnitude of these effects: even 
for K = \ J\ these are in the percent range, and for the 
physical parameters of cuprate systems, if ~ J/5, they 
are quantitatively irrelevant. 



B. Projected model 

In Sec. IV we have shown that circulating-current 
states are significantly more competitive in projected 
models for strongly correlated electrons than in HFA. On 
first inspection one may also expect ring-exchange inter- 
actions to be considerably more relevant in a projected 
Hamiltonian: the appropriate Gutzwiller prefactor for 
the K term is gK = gj, which approaches 16 at half- 
filling. To examine this possibility, in Fig. 10 we show 
for an intermediate doping the projected analog of Fig. 8, 
following the bond current and ring-exchange energy in a 
fully projected CF state as a function of K. Both quanti- 
ties show the same essentially linear trend which follows 



from suppression of the CF phase in the presence of a 
positive ring-exchange interaction. With regard to the 
doping-dependence of Icf and Ek.cf, as in the unpro- 
jected case the current is dictated by the K = behavior 
of the system [Fig. 6(a)], with a small reduction in criti- 
cal doping for positive K. The ring-exchange energy for 
the projected model with K = J/5 is shown in Fig. 9(b), 
and exhibits a steady decline which may be traced to the 
combination of parameters [Eqs. (12, 13)] as the bond 
phase decreases from 7r/4 to [Fig. 6(b)]. 

Returning to Fig. 10, in fact the ring-exchange effect 
on the current in the projected CF phase is significantly 
larger than in the unprojected one (10% compared to 2% 
at K = J for X = 0.083). This result holds for all dop- 
ings, and extends also to the relative magnitude of the 
energy contribution Ex- However, it is clear that changes 
in both energy and currc;nt remain extremely small: this 
is not surprising given that the energetic contribution of 
the K term is multiplied by the fourth power of bond 
expectation values Sij which are suppressed towards zero 
at half-filling for full projection. In this context we note 
that ordered magnetic states suffer no such suppression 
of the order parameter on projection, and thus from the 
large value of gK a Gutzwiller-projected model near half- 
filling might be expected to favor the formation of mag- 
netic phases stabilized by large values of K. 



C. Phcises in the large-K regime 

We thus conclude this section by commenting briefly 

on the effects of a very strong ring-exchange interac- 
tion. In spin systems, as realized at half-filling for 
large U, the ring-exchange term is generally considered 
to favor a maximization of the solid angle of the four 
spins on a plaquette, and thus to yield non-coUinear 
magnetic phases. In the large-if limit one expects the 
(static^^ or dynamic^^'^^) vector-chirality phase depicted 
in Fig. 11(a). This type of spin configuration has a low 
bond kinetic energy but a finite spin fiux in the sense fa- 
miliar from double-exchange models,^^ that a propagat- 
ing electron picks up a path-dependent phase. In practice 
the bond current and kinetic energy vanish at half-filling 
in this static limit for large-?/ systems. In the presence of 
holes a strong ring-exchange interaction is also expected 
(from the magnitude of the expectation values involved) 
to promote magnetic rather than bond order. 

However, in our analysis of the large- ii' limit for an un- 
projected model with U/t » 1, both at and away from 
half-filling, we have found a different ground state, shown 
in Fig. 11(b). This spin configuration ensures a "three- 
up, one-down" state on every plaquette, which gives the 
same ring-exchange energy per plaquette [see Eq. (4)] as 
the vector-chirality configuration. It is not unique, and 
need not conserve a total spin component 5^°* = 0. The 
extra stability of the "three-up, one-down" phase may be 
ascribed primarily to the fact that, in approaching the 
limit of an effective spin model, residual dynamical pro- 
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FIG. 11: (a) Schematic representation of a static spin con- 
figuration with vector-chirality order, (b) A "threo-up, one- 
down" configuration, as obtained for x = both in an unpro- 
jected model using JJ = lOt, J = W = O.Olt, and K = 25t, 
and in a fully projected model with J = 4V = O.Olt and 
K = Wt. 

cesses act to favor coUinear spin states. However, quan- 
tum fluctuations not contained at the Hartree-Fock level 
could also play an important role in the stability of the 
phases in Fig. 11. 

In the projected model obtained by consistent deriva- 
tion from the large-f/ Hubbard Hamiltonian, we have 
been unable to find these states because of continued 
competition from supcrexc;liange terms, which are always 
significant. In this case the mutual frustration of the two 
competing interactions raises the possibility of further 
exotic spin configurations, such as the scalar-chirality 
phase, or different dimerized structures breaking lat- 
tice translation symmetry. ^^"^^ However, with an arbi- 
trary choice of large K unrelated to the consistency con- 
siderations of Sec. IV, the same "three-up, one-down" 
phase emerges directly [Fig. 11(b)]. 

VI. SUMMARY AND CONCLUSIONS 

In this contribution we have considered a microscopic 
description of circulating charge- and spin-current states 
in a planar model appropriate for cuprate systems, and 
the qualitative influence upon these phases of cyclic, 
four-spin interactions. The CF and SF phases were 
established on a small cluster by iterative solution for 
the electronic states within the HFA for the general- 
ized Hubbard model. Although the real-space Hartree- 
Fock approach does not in general favor homogeneous 
and dynamically stabilized states, prefering localized and 
charge-inhomogeneous structures, meaningful qualitative 
comparisons remain possible. 

The robust qualitative conclusions concerning fiux 
phases are the following. CF and SF states are favored 
for small on-site repulsion U, and are stabilized primarily 
by next-neighbor supcrcxchange (J) and Coulomb {V) 
contributions. They are most favorable close to half- 
filling, becoming unstable to non-ordered metallic phases 
at higher hole doping (x > 0.15 in HFA), a result which 
refiects qualitatively the onset of the Fermi-liquid phase 
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beyond optimal doping in the cuprates. The bond cur- 
rent / increases with interaction strength, tending in the 
limit |J|/t, 4V/t 1 towards a 7r-fiux phase in the CF 
state, and in the SF state to an imaginary bond order pa- 
rameter preceding a ferromagnetic instability. Very few 
inhomogeneous current states are found in this type of 
analysis, and coexisting states are dominated by local- 
ized charge and magnetic configurations. 

The inclusion of strong electronic correlations by pro- 
jection of the hopping terms may yield solutions qual- 
itatively different from those obtained for independent 
electrons in a weak-coupling treatment (unprojected hop- 
ping). Projected hopping is therefore essential in making 
contact with cuprate-relevant parameters. Inclusion of 
projection effects within the Gutzwiller approximation 
shows that the CF phase does become competitive for 
realistic cuprate interactions. The CF phase in the pro- 
jected model has a clearly defined regime of stability at 
low and intermediate doping, in close correspondence to 
the underdoped and optimally doped regimes in cuprates. 
While projection produces very small kinetic terms near 
half-filling, it also pushes the general CF state towards 
a TT-flux phase. Next-neighbor hopping contributions (t') 
of the magnitude and sign relevant for cuprate systems 
appear to reinforce the projected CF phase. By contrast, 
the SF phase is excluded in a projected model, demon- 
strating that this type of flux phase is not compatible 



with strong electron correlations. 

Introduction of a ring-exchange interaction {K > 0) is 
detrimental to the stability of CF phases in both unpro- 
jected and projected models. By contrast, the SF phase 
shows an enhanced current and minimal energy penalty 
in the unprojected framework. Ring-exchange eff'ects 
vary little with doping until the spin and bond order pa- 
rameters are strongly reduced at high filling. The param- 
eter values in cuprates are such that the K term is largely 
irrelevant for qualitative phenomena at the Hartree-Fock 
level, and its effect remains weak even with the strong 
enhancement provided by the Gutzwiller approximation. 
However, we have shown that an exception to this state- 
ment arises for ordered magnetic states at very low fill- 
ing, where the ring-exchange interaction might become 
qualitatively significant in a projected model. 
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